Topological Phase Transitions in the Golden String-Net Model 
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We examine the zero-temperature phase diagram of the two-dimensional Levin- Wen string-net 
model with Fibonacci anyons in the presence of competing interactions. Combining high-order 
series expansions around three exactly solvable points and exact diagonalizations, we find that the 
non-Abelian doubled Fibonacci topological phase is separated from two trivial phases by different 
second-order quantum critical points, the positions of which are computed accurately. The trivial 
phases are separated by a first-order transition occurring at a fourth exactly solvable point where the 
ground-state manifold is infinitely-many degenerate. The evaluation of critical exponents suggests 
unusual universality classes. 
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Quantum phases of matter are often well described by 
local order parameters and Landau-Ginzburg symmetry- 
breaking theory is an efficient tool to analyze transi- 
tions between these phases. However, in the late 1980s, 
a new class of phases that cannot be understood in 
terms of local symmetries has emerged in the context of 
high-temperature superconductivity [1-3]. These phases, 
dubbed topological because of their sensitivity to the sys- 
tem topology, have stimulated many studies in different 
domains (see Ref. 4 for a recent review) . One of the most 
intriguing properties of topologically-ordered systems is 
that they are robust against local (not too strong) pertur- 
bations [5, 6]. This stability makes them especially ap- 
pealing for quantum computation [ ] as well as good can- 
didates for quantum memories [8]. Several experiments 
have been proposed to realize the so-called topologically- 
protected qubits [9]. In this perspective, a theoretical 
characterization of the robustness of topological phases 
under strong perturbations as well as the nature of the 
phase transitions signaling their breakdown is undoubt- 
edly an important issue. Thanks to recently proposed 
exactly solvable lattice models realizing various topolog- 
ical phases of matter [5, 10, 11], this program has been 
undertaken for several models [12-21]. 

The main purpose of this Letter is to go one step 
beyond, by studying the phase diagram of a paradig- 
matic 2D non-Abelian model. We consider the Levin- 
Wen model [11] on the honeycomb lattice with Fibonacci 
anyons (the golden string-net) in the presence of the same 
perturbation as the one introduced in Ref. 16. We de- 
termine the extension of the doubled Fibonacci (DFib) 
topological phase and show that it is separated from two 
other topologically trivial phases via second-order tran- 
sitions that are analyzed in details. 

Hilbert space — Following the Levin- Wen construction, 
we consider a honeycomb lattice with anyonic degrees of 
freedom living on its edges. In the Fibonacci string-net 



model, these local (microscopic) degrees of freedom can 
be in two different states |0) or |1). The Hilbert space H 
is restricted to states that satisfy the so-called branching 
rules stemming from the non-Abelian fusion rules 

x a = a x = a for a € {0, 1}, (1) 
1x1 = 0+1. (2) 

At each vertex of the honeycomb lattice, the fusion rules 
must not be violated, i.e., if one edge is in state |1), then 
at least one of the two other edges must be in the same 
state. For an arbitrary trivalent graph with N v vertices, 
the dimension of the Hilbert space is then given by 

dim'H = (l + ^ 2 )^ 1 +(l + <^ 2 )^\ (3) 

where ip = is the golden ratio (see, e.g., Ref. 22). 

Model — We study the following Hamiltonian [16] 

H = ~Jp 2J ^*(P).0 _ Jc 2J <^(e),0- (4) 
p e 

The first term is the string-net Hamiltonian introduced 
by Levin and Wen [11]. It involves the projector 5$( p ) ; o 
onto states with no flux $(p) through plaquette p [11, 16]. 
The second term is diagonal in the basis introduced above 
since <5;( e ),o is the projector onto state |0) on edge e. 

To help the reader to grasp the physical content of 
this Hamiltonian, let us mention what happens if one 
replaces Eq. (2) by the simpler Abelian Z2 fusion rule 
1x1 = 0. The model then becomes Kitaev's toric code 
[ ] on the honeycomb lattice, restricted to the charge-free 
sector (because of the branching rules), in the presence 
of a magnetic field J c in the x direction. Indeed, one 
can write the Hamiltonian in terms of Pauli matrices, 
with 8i<j, )fl = (1 + rUp CT e)/2 and 5 l(e)fi = (l + of)/2. 
Identifying plaquette fluxes with spin 1/2 variables as 
done in Refs. 12 and 13 for the square lattice, the Hamil- 
tonian can further be mapped onto the transverse-field 
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Ising model on the triangular lattice (of plaquettes) , with 
coupling J G and transverse-field J p . 

Limiting cases — For convenience, let us set J p = cos 9 
and J c = sin 9. To our knowledge, contrary to the lad- 
der geometry studied in Ref. 16, the Hamiltonian H is 
exactly solvable only at the four points for which 9 is 
a multiple of ir/2. In the following, we discuss the low- 
energy spectrum of H as well as the corresponding phases 
around these special points and give some arguments in 
favor of transitions between them. 

• = 0: For 6 = 0, the model reduces to the golden 
string-net model. Ground states |g) are flux- free states 
satisfying 5$(p),o|g) = |g) for all p and thus have an en- 
ergy per plaquette eo = Eq/N p = — 1. Their degeneracy 
depends on the system topology which is the most salient 
property of a topological phase. For the Fibonacci string- 
net model on any trivalent graph, the ground state is 
unique on a sphere [23] whereas it is four-fold degenerate 
on a torus. Interestingly, one can also compute the degen- 
eracy of the k th excited states (with energy = E + k) 




where we introduced the number of plaquettes N p and 
the famous Fibonacci sequence defined for any integer n 
by F n+1 = F n + F„_i, with F-i = 1 and F = 0. The 
integers (p, q, r) depend on the surface considered. For 
instance, one has (p, q, r) = (1, 0, 0) on a sphere whereas 
(p,q,r) = (4,1,4) on a torus. Equation (5) shows that 
for the Fibonacci theory, an odd number of excitations 
can exist on a compact surface contrary to the charge- 
free toric code discussed above where fluxes are always 
created and annihilated by pairs. Note that the binomial 
coefficient simply arises from the different ways to choose 
k plaquettes carrying the flux excitations among N p . 

Products of two iVs stem from the fact that the "emer- 
gent" flux excitations are not the microscopic Fibonacci 
anyons but are achiral combinations of two Fibonacci 
anyons [ ]. The non-Abelian topological phase in the 
vicinity of 9 = is described by a DFib theory [11, 25- 
29]. Excitations have a trivial topological spin [11] and 
can fuse to the vacuum (also called trivial particle) [30]. 
As such, they can also be considered as bosons [31] and 
hence condense. 

• 9 = tt : For J c = and J p < 0, the low-energy spec- 
trum is very different. Indeed, in this case, the ground- 
state manifold is 2?7v p -fold degenerate and spanned by 
all states \g) satisfying £>$( p ),o|g) = for all p. As dis- 
cussed above, this degeneracy depends on the topology 
through its indices (p, q, r) so that one may still consider 
the system as topologically ordered. 

Owing to this huge degeneracy, we have not been able 
to analyze the vicinity of this point. However, numerical 
results obtained by exact diagonalizations clearly show 
that (i) the degeneracy is lifted as soon as the coupling 
J c 7^ and (it) the ground state for 9 = tt^ is unique and 
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FIG. 1. (Color online). Ground-state energy per plaquette 
eo = Eo/Np as a function of 0. ED results (black line) for 
N p = \/T3 x vT3 plaquettes (see inset) are in excellent agree- 
ment with typical Pade approximants (white lines) computed 
from high-order series expansions [32]. 

adiabatically connected to the polarized ground states 
found at 9 = tt/2 and 37r/2 respectively (see discussion 
below). This result is in stark contrast with the scenario 
described in Ref. 1 6 on the ladder where a gapless phase 
is observed for 9 € [7r,37r/2]. In addition, as can be 
seen in Fig. 3 (central panel), we found a jump in deeo 
at 9 = tt (for all system sizes) indicating that the two 
gapped phases (9 = tt + and 9 = n~) are separated by 
a first-order phase transition. This means that there is 
no topological phase beyond the point 9 = tt. In some 
sense, one may say that the topological order present at 
this point is extremely fragile. 

• 9 = tt/2 : For J p = 0, the Hamiltonian H is diagonal 
in the canonical basis of states satisfying the branching 
rules. For J c > 0, the ground state is unique whatever 
the topology and corresponds to the fully polarized state 
where all edges are in the state |0) (with eigenenergy 
Eq = —N e where N e is the total number of edges). First 
excited states are obtained by flipping six links around 
one hexagon. They behave as trivial hard-core bosons 
that become dynamical when the coupling J p is switched 
on. Thus, near 9 = tt/2, the system is gapped but not 
topologically ordered, making the occurrence of a phase 
transition in the interval [0, 7r/2] compulsory. 

• 9 = 3tt/2 : For J < and J p = 0, the Hamil- 
tonian is also diagonal and the unique ground state is 
the fully polarized state where all edges are in the states 
|1) (eo = 0). Note that such a state would be forbidden 
by the Abelian Z2 fusion rules. First excited states are 
obtained from the ground state by flipping a single link. 
As previously, these localized excitations are trivial hard- 
core bosons that become mobile when J p 7^ so that one 
expects a phase transition in the interval [37r/2, 27r]. 

Phase diagram — To determine the zero-temperature 
phase diagram, we combined two different approaches. 
First, we performed high-order series expansions in the 
thermodynamical limit by means of several methods [33- 
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FIG. 2. (Color online). First four excitation energies obtained 
from ED results (black lines) for N p = 3x3 plaquettes com- 
pared with the low-energy gap computed from high-order bare 
series expansions (white lines) [32]. The topological degener- 
acy splitting is clearly observed in the vicinity of the critical 
points. For symmetry reasons, this splitting is only partial 
for the system considered here (see inset). 

.], around the exactly solvable points 9 = 0, ir/2, and 
37r/2 described above. This yields the ground-state en- 
ergy per plaquette eo as well as the quasiparticle disper- 
sion from which the low-energy gap A can be extracted. 
Lengthy expressions of these series expansions can be 
found in Ref. 32. This method allows one to accurately 
compute the critical couplings for which the gap van- 
ishes. These points are associated to second-order tran- 
sitions but might not be relevant if first-order transitions 
are present (see Refs. 36 and 37 for details about this 
issue in a similar context). Second, we perform exact di- 
agonalizations (ED) for lattices with periodic boundary 
conditions. As can be seen in Figs. 1 and 2, series ex- 
pansions and ED data are in very good agreement except 
in the vicinity of the transition points where finite-order 
and finite-size effects are important. 

Combining these two methods (ED and series expan- 
sions) we found that the DFib topological phase near 
9 = ranges from 0§ ~ -0.63 (= 5.65) to 0\ ~ 0.255. 
As we shall now argue, we associate these two points 
to second-order transitions. The first piece of evidence 
pleading in favor of such a scenario follows from ED and 
is the behavior of <9 2 eo that clearly decreases with the 
system size near these points (see left and right panels 
in Fig. 3). In addition, the position of the low-energy 
gap minimum as well as the topological degeneracy split- 
ting shown in Fig. 2 lie in the same region as the po- 
sition of the minimum of c^erj. Let us also note that 
we did not find any relevant level crossing in the exci- 
tation spectrum that could lead to a first-order transi- 
tion. The second argument comes from the high-order 
perturbation theory. As can be seen in Ref. 32, the se- 
ries behave very differently for positive and negative J e . 
So, let us first discuss the most favorable case J c > 



(J p > 0). A close inspection of the series expansion near 
9 = and 9 = ir/2 reveals three important features: (?) 
the sign of each coefficient is the same in all series ; (it) 
series of the ground-state energy intersect in two order- 
dependent points - the possible merging of these points, 
in the infinite-order limit, being a signal of a second-order 
transition since series then have to be tangential at the 
critical point ; (Hi) series of the gap intersect in a unique 
(still order-dependent) point, see, e.g., left part of Fig. 2. 
We emphasize that the value of the gap at this crossing 
point decreases when the order increases and eventually 
vanishes in the infinite-order limit. In Fig. 4, we plot- 
ted the position of these crossing points as a function of 
the (inverse) order as well as the position of the mini- 
mum of the low-energy gap and of (9 2 eo as a function of 
Np 1 computed from ED results. As can be seen, all data 
seem to converge to the same point 9\ E [0.255,0.256] in 
the infinite-order (size) limit, providing a smoking-gun 
evidence of a second-order transition. 

Unfortunately, the case J c < (J p > 0) is more in- 
volved for three reasons. First, series expansions of eo 
and A in this region have alternate signs so that the pre- 
vious criteria based on crossing points cannot be used. 
Second, contrary to the case J c > 0, the momentum min- 
imizing the dispersion of the low-energy quasiparticles is 
not at the T point and only belongs to the reciprocal lat- 
tice of 3 p x 3 q systems, (p, q) € N 2 . The only system with 
such characteristics considered in this study is the 3x3 
lattice so that one cannot perform any reliable analysis 
from ED data near 9\. Third, because of the nature of 
the low-energy states, maximum orders reachable around 
9 = 3ir/2 are smaller than around 9 = ir/2 [32]. Never- 
theless, using the gap series expansion around 9 = 0, 
it is possible to perform (DlogPade) resummations that 
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FIG. 3. (Color online) ED results for N p = 2 x 2 (dotted line), 
3x3 (dashed line), and \/l3 x \/T3 (solid line) plaquettes. 
Left and right panels : c*f eo decreases with the system size in- 
dicating second-order transitions at 91 and Q\ ■ Central panel : 
dgeo displays a clear jump at — n indicating a first-order 
transition. Dips indicated by arrows in the DFib topological 
phase are due to (irrelevant and avoided) level crossings be- 
tween the four lowest-energy levels that become degenerate 
in the thermodynamical limit on a torus. 
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FIG. 4. (Color online) Position of different quantities as a 
function of either the inverse order n _1 of the corresponding 
series expansion [crossing points of eo (dots) , crossing points 
of the gap (triangles)] or Np 1 [minimum of the gap (squares), 
minimum of <9|eo (diamonds)]. Lines are power-law fits ob- 
tained by choosing 81 such that it maximizes the correlations 
for the crossing points of eo (solid lines) that are the most ac- 
curate results in this work. These crossing points have been 
obtained from order n series around 8 = and order 2n — 1 
series around 8 = tt/2. For the crossing points of the gap, we 
used order n series around both limits. 



lead to a position of the critical point 9% e [5.61,5.62]. 
Note that the same methods for J c > would lead to 
9\ e [0.259,0.261]. This value differs from the one pro- 
posed above only by a few percents and it overestimates 
the extension of the topological DFib phase. Conse- 
quently, as we have no alternative approach, we roughly 
estimate that 9\ lies in the range [5.6,5.7]. 

Critical exponents — The obvious question that arises 
next concerns universality classes associated to the tran- 
sition points. In the absence of a local order parame- 
ter, the only meaningful critical exponents for topological 
phase transitions are those associated to the spectrum. 
Let us remind that for a second-order transition, the gap 
vanishes, at large linear system size L and at the critical 
point, as A ~ L~ z where z is the dynamical exponent. In 
the thermodynamical limit, one further has A ~ \9— 9 C \ ZU 
and dje ~ |0-0 c |- a . 

As already explained, ED are only useful quantita- 
tively around 9\ (although restricted to systems of small 
sizes and different shapes). Using the gap data, we find 
z ~ 1.2. A finite-size analysis of eo yields a surprisingly 
good data collapse for the 3x3 and vT3 x \/T3 sys- 
tems, with 9\ ~ 0.255, z 1, v ~ 0.4 and a ~ 0.8 [24]. 
The rather large value of a might be responsible for the 
quality of the data collapse. We emphasize that these 
values are compatible with the previous estimate of 9\ as 
well as with the hyperscaling relation 2 — a = 2/(2 4- z). 
The above exponents are furthermore in agreement with 
DlogPade resummations of the series expansion around 
9 = which yield zv £ [0.4,0.44]. We discard results 
from resummations around 9 — ir/2 because these re- 
sults are not as well converged, since for a given order, 



clusters needed to compute series around this point are 
twice as small as those needed around 9 = 0. Note that 
as usual, extracting a from series of eo does not give any 
conclusive result. Concerning the critical behavior at 9%, 
we only use DlogPade resummation around 9 = and we 
obtain a gap exponent zv € [0.56,0.58]. 

Outlook — It is difficult to provide some error bars 
concerning these values. We benchmarked our series ex- 
pansion analysis on the Fibonacci ladder (for which ex- 
ponents are known exactly [16]) and the 2D Z2 string-net 
model (having either Ising or XY transitions depending 
on the sign of J c ). The results we obtained [24] lead 
us to conclude that critical exponents are to be consid- 
ered with a precision of about 10%. As a conclusion, we 
found two different second-order transitions with univer- 
sality classes that, to the best of our knowledge, are as 
yet unknown in the context of topological phase transi- 
tions. Let us mention that a critical DFib wavefunction 
have been proposed [38] but its relevance for the present 
problem require further studies [ ] . 

To gain more understanding about these transitions, 
different approaches could be used, e.g., variational meth- 
ods or Monte-Carlo simulations, although a naive im- 
plementation of the latter should suffer from the sign 
problem. It would also be worth studying similar models 
with a DFib phase [38, 39] as well as different topological 
phases. Finally, another important issue concerns tran- 
sitions between two distinct topological phases [3 ]. We 
hope that the present work will stimulate such investiga- 
tions. 
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SUPPLEMENTAL MATERIAL 

In the following, we give the series expansions in the different phases for the ground-state energy per plaquette 
eo and the quasiparticle gaps A , for positive and negative signs of the dimensionless parameter t — J c /J p = tan# 
respectively. For the sake of clarity, we give below the numerical values of the coefficients with 16 digits. 



Expansions in the vicinity 6 — 

The ground-state energy per plaquette e near 9 — (J p = 1, J e = 0) has been obtained up to order 10 using 
operator perturbation theory [34], whereas quasiparticle gaps A* were obtained up to order 9 using perturbative 
continuous unitary transformations. [35] 



eo/Jp = - 1 - 0.8291796067500631 i - 0.3 t 2 - 0.2329179606750063 t 3 - 0.3758359213500126 t 4 

- 0.6934622369921362 1 5 - 1.517757831138397 1 6 - 3.615896887905089 t 7 - 9.257482947753094 1 8 

- 24.89646210135949 t 9 - 68.81054927770635 1 10 , (6) 

A+/J p = 1 - 1.658359213500126 1 - 2.029179606750063 t 2 - 3.107113095525145 t 3 - 8.042597266963313 t 4 

- 19. 16249885423558 1 5 - 58.31720409052607 1 6 - 164.4421257647495 t 7 - 528.5318111014412 1 & 

- 1620.571984517245 t 9 , (7) 



A"/J p = 1 + 0.8291796067500631 1 + 0.1145898033750315 1 2 + 0.5110332556124590 t 3 + 0.4044760408194119 t 4 
+ 0.9554382719956335 1 5 + 1.784752477740017 1 6 + 4.523961920423115 1 7 + 11.17294663306187 1 8 
+ 30.67329262671686 1 9 . (8) 



Expansions in the vicinity 8 = n/2 

The ground-state energy per plaquette e near 9 = 7r/2 ( J c = 1, J p = 0) has been obtained up to order 19 using a 
partitioning technique provided by Lowdin [33], whereas quasiparticle gaps A ± were obtained up to order 11 using 
operator perturbation theory [34] on appropriate periodic clusters. 

e /J c = - 3 - 2.763932022500210 • 10" 1 r 1 - 3.333333333333333 • 10~ 2 r 2 - 2.484519974999766 • 10~ 3 1~ 3 



- 3 - 2.763932022500210 


10" 




1 - 3.333333333333333 • 10~ 2 


r 2 - 2.484519974999766 


• 10 


- 1.473090114646592 


10" 


4 r 


4 _ 


1.762516450320833 • 10~ 5 <~ 5 


- 3.114829150546602 • 10 


" 6 <- 


- 4.974954832303385 


10" 


7 r 


7 _ 


8.712942025753695 • 10~ 8 i~ 8 


- 1.680470831303724 • 10 


" 8 <- 


- 3.252798060742993 


10" 


9 r 


10 _ 


- 6.452247707667803 • lO -10 *" 


11 - 1.328942917494399- 


10 - 


- 2.786531244027440 


10" 


11 1 


-13 


- 5.923931315552463 • 10~ 12 1 


~ 14 - 1.280290681522505 


• 10" 


- 2.803432623420322 


10" 


13 1 


-16 


- 6.202995686227467 • 10~ 14 t 


" 1T - 1.386248233596245 


• 10" 


- 3.125540962122164 


10" 


15 1 


-19 









10 1- 12 

-12 f -15 
-14 t -18 

(9) 

A ± /J c = 6 - 4.472135954999579 • 10" 1 r 1 - 5.916925468334595 • 10~ 2 r 2 - 1.532368203449658 • 10~ 2 r 3 

- 1.826128841179490 • 10~ 3 i~ 4 - 2.141218566177655 ■ 10~ 4 i" 5 - 7.057273308076495 • 10~ 5 ^ 6 

- 7.713120260517609 • 10~ 6 ^ 7 - 1.822619200222926 ■ 10~ 6 i" 8 - 1.739394034185955 • 10~ 7 t~ 9 

- 1.068512192171555 • 10~ 7 r 10 - 1.583732801424164 • 10~ 8 t~ n . (10) 
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Expansions in the vicinity 9 — 3n/2 

The ground-state energy per plaquette eo near 9 = 3n/2 (J e = — 1, J p = 0) has been obtained up to order 9 using 
partitioning techniques [33], whereas quasiparticle gaps A were obtained up to order 6 using perturbative continuous 
unitary transformations. [35] 

eo/ (~J e ) = 3.013155617496425 • 1CT 1 t' 1 - 1.132044933254820 • 10" 1 r 2 + 2.807797460719963 • 10~ 2 r 3 

- 4.507809490972001 • 10~ 3 tr 4 - 3.037988446794379 • 10" 3 i~ 5 + 4.596752532754356 • 10~ 3 1~ 6 

- 1.633669374413878 • 10 -3 1~ 7 - 2.288717115668441 • 10~ 3 i~ 8 + 3.812642130493073 • 10~ 3 1~ 9 , (11) 

A+/ (- J e ) = 1 - 7.331262919989905 • 10" 1 t~ l + 2.627515502196252 • 10" 1 r 2 - 1.016556730679853 • lO^ 1 t' 3 

+ 1.966652521926989 • 10~ 2 r A - 7.070423645089451 ■ 10~ 2 1~ 5 + 1.098516898950685 • 10" 1 1~ 6 , (12) 

A"/ (- J c ) = 1 + 2.668737080010095 • 10" 1 1~ 1 - 2.048165955279984 • 10" 1 r 2 - 7.571315893653220 • 10~ 2 t~ 3 

+ 3.964262957878962 • 10~ 2 t~ A + 7.257956045336907 • 10~ 2 i~ 5 - 6.691149658050659 • 10~ 2 r 6 . (13) 



